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DESCRIPTION 

CROSS REFERENCE TO RELATED APPLICATIONS 

This application is a continuation of International 
Application Number PCT/GB/03073 , filed 3 July 2003, which is 
incorporated by reference herein. 

TECHNICAL FIELD 

The invention relates to contact sensitive devices. 
BACKGROUND ART 

Visual displays often include some form of touch 
sensitive screen. This is becoming more common with the 
emergence of the next generation of portable multimedia 
devices such as palm top computers. The most established 
technology using waves to detect contact is Surface Acoustic 
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Wave (SAW) , which generates high frequency waves on the 
surface of a glass screen, and their attenuation by the 
contact of a finger is used to detect the touch location. This 
technique is "time-of -flight" , where the time for the 
5 disturbance to reach one or more sensors is used to detect the 
location. Such an approach is possible when the medium behaves 
in a non- dispersive manner i.e. the velocity of the waves does 
not vary significantly over the frequency range of interest. 

In contrast in WO01/48684 to the present applicant, a 

10 contact sensitive device and method of using the same are 
proposed. The device comprises a member capable of supporting 
bending wave vibration and a sensor mounted on the member for 
measuring bending wave vibration in the member and for 
transmitting a signal to a processor whereby information 

15 relating to a contact made on a surface on the member is 
calculated from the change in bending wave vibration in the 
member created by the contact. 

By bending wave vibration it is meant an excitation, for 
example by the contact, which imparts some out of plane 

20 displacement to the member. Many materials bend, some with 
pure bending with a perfect square root dispersion relation 
and some with a mixture of pure and shear bending. The 
dispersion relation describes the dependence of the in-plane 
velocity of the waves on the frequency of the waves. 

25 Two types of contact sensitive device are proposed, 

namely a passive sensor in which bending wave vibration in the 
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member is only excited by the contact and an active sensor in 
which the contact sensitive device further comprises an 
emitting transducer for exciting bending wave vibration in the 
member to probe for information relating to the contact. In 
5 the active sensor, information relating to the contact is 
calculated by comparing the response of waves generated by the 
emitting transducer in the absence of a contact to the 
response caused by the mechanical constraint of the presence 
of a contact . 

10 Bending waves provide advantages, such as increased 

robustness and reduced sensitivity to surface scratches, etc. 

However, bending waves are dispersive i.e. the bending wave 
-velocity, and hence the "time of flight", is dependent on 

frequency. In general, an impulse contains a broad range of 
15 component frequencies and thus if the impulse travels a short 

distance, high frequency components will arrive first. This 

effect must be corrected . 

In WO01/48684, a correction to convert the measured 

bending wave signal to a propagation signal from a non- 
20 dispersive wave source may be applied so that techniques used 

in the fields of radar and sonar may be applied to detect the 

location of the contact. The application of the correction is 

illustrated in Figures la to Id. 

Figure la shows an impulse in an ideal medium with a 
2 5 square root dispersion relation and demonstrates that a 

dispersive medium does not preserve the waveshape of an 
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impulse. The outgoing wave (60) is evident at time t=0 and the 
echo signal (62) is spread out over time, which makes a 
determination of an exact contact position problematic. 

In a non-dispersive medium such as air, a periodic 
5 variation of the frequency response is characteristic of a 
reflection, and is often referred to as comb filtering. 
Physically, the periodic variation in the frequency response 
derives from the number of wavelengths that fit between the 
source and the reflector. As the frequency is increased and 

10 the number of wavelengths fitting in this space increases, the 
interference of the reflected wave with the outgoing wave 
oscillates between constructive and destructive. 

Calculating the Fourier transform of the dispersive 
impulse response of Figure la produces the frequency response 

15 shown in Figure lb. The frequency response is non-periodic 
and the periodic variation with wavelength translates to a 
variation in frequency that gets slower with increasing 
frequency. This is a consequence of the square root 

dispersion in which the wavelength is proportional to the 

2 0 square root of the inverse of frequency. The effect of the 
panel on the frequency response is therefore to stretch the 
response as a function of frequency according to the panel 
dispersion. Consequently, a correction for the panel 
dispersion may be applied by applying the inverse stretch in 

2 5 the frequency domain, thus restoring the periodicity present 
in the non-dispersive case. 
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By warping the frequency axis with the inverse of the 
panel dispersion, Figure lb may be transformed into the 
frequency response for the non-dispersive case (Figure lc) in 
which the frequency of excitation is proportional to the 
5 inverse of the wavelength. This simple relationship translates 
the periodic variation with decreasing wavelength to a 
periodic variation with increasing frequency as shown in 
Figure lc . 

Applying the inverse Fast Fourier Transform (fft) to the 
10 trace of Figure lc produces an impulse response shown in 
Figure Id which is corrected for dispersion and where the 
clear reflection is restored. As is shown in Figure Id any 
particular waveshape of an impulse is preserved in time since 
the waves travelling in a non-dispersive medium have a 
15 constant velocity of travel, independent of their frequency. 
Accordingly, the task of echo location is relatively straight 
forward. The outgoing wave (50) is evident at time t = 0, 
together with a clear reflection (52) at 4ms. The reflection 
(52) has a magnitude which is approximately one-quarter of the 
2 0 magnitude of the outgoing wave (50) . 

The procedure described is not applicable if the impulse 
has occurred at an unknown time t 0 and the distance x from the 
response to an initial impulse may only be calculated if the 
impulse occurs at t 0 = 0 . 

25 It is an object of the present invention to provide an 

alternative contact sensitive device which uses bending wave 
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vibration for extracting information relating to the contact. 
DISCLOSURE OF INVENTION 

According to one aspect of the invention, there is 
provided a contact sensitive device comprising a member 
5 capable of supporting bending waves, a first sensor mounted on 
the member for measuring bending wave vibration in the member, 
the first sensor determining a first measured bending wave 
signal and a processor which calculates information relating 
to a contact on the member from the measured bending wave 

10 signal, the processor applying a correction based on the 
dispersion relation of the material of the member supporting 
the bending waves, wherein the device comprises a second 
sensor to determine a second measured bending wave signal 
which is measured simultaneously with the first measured 

15 bending wave signal and the processor calculates a dispersion 
corrected function of the two measured bending wave signals 
which is selected from the group consisting of a dispersion 
corrected correlation function, a dispersion corrected 
convolution function, a dispersion corrected coherence 

2 0 function and other phase equivalent functions to determine 
information relating to the contact . 

According to a second aspect of the invention, there is 
provided a method of determining information relating to a 
contact on a contact sensitive device comprising the steps of 

25 providing a member capable of supporting bending waves and a 
first sensor mounted on the member for measuring bending wave 
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vibration in the member, determining, using the sensor, a 
first measured bending wave signal, providing a second sensor 
mounted on the member to determine a second measured bending 
wave signal, measuring the second measured bending wave signal 
5 simultaneously with the first measured bending wave signal, 
calculating a dispersion corrected function of the two 
measured bending wave signals which is selected from the group 
consisting of a dispersion corrected correlation function, a 
dispersion corrected convolution function, a dispersion 
10 corrected coherence function and other phase equivalent 
functions and processing the measured bending wave signals to 
calculate information relating to the contact by applying the 
dispersion corrected function. 

The following features may be applied to both the device 
15 and the method with the processor being adapted to provide 
many of the calculations or processing steps of the method. 

The dispersion corrected function may be calculated as 
follows : 

calculate W x (a>) and W 2 (a))* which are the Fourier 
20 transformation and complex conjugate Fourier transformation of 
the two measured bending wave signals W x {t) and W 2 {t) ; t 
represents time co is 2nt where f is frequency. 

calculate a first intermediate function W x {co)W 2 \co) ; 
calculate a second intermediate function M(o) which is a 
25 function of W x {(oW 2 {co) ; 

Atty. Dkt. No. 081-0002-1 
081-0002-1 final patent app.doc 



8 

apply a frequency stretching operation f(co), as 
described above in relation to WOO 1/48684 , to M(co) to give the 
dispersion corrected correlation function: 



G(t) = — rM[f(a))]exp(iG)t)dco . 

9 rr 



In 

5 The intermediate function M(co) may simply be W x {&}W*(cS) 

which gives a standard dispersion corrected correlation 
function. Alternatively, M(ca) may be a function which modifies 
the amplitude but not the phase of W x (a>)W*(a>) to give a phase 
equivalent function to the standard dispersion corrected 

10 correlation function. Since the phase equivalent function and 
the standard dispersion corrected correlation function have 
the same phase properties, they have a maximum at the same 
position. The phase information in the measured bending wave 
signals may be used to acquire information about the contact 

15 in particular the location thereof. The location may be 
calculated from the time at which the maximum in the functions 
occurs . 

M(co) may be selected from the following functions which 
all yield phase equivalent functions to the standard 
2 0 dispersion corrected correlation function: 

W x (cq)W;(g>) 



W x (co)W 2 \co)\ 



a) M(co) = 



Thus M(to) may normalise the amplitudes of W } (co)W*(a>) to 
unity to yield a normalised dispersion corrected correlation 
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function otherwise known as a dispersion corrected coherence 
function. 



b) MW . twKw 



Thus M(co) may act on the amplitudes of W x (a>)W* (a>) to 
5 yield a dispersion corrected correlation function with a 
modified peak shape. 

c) M(co) = ^(^)^(^)^[#j(^)^(^)|J where (p{x) is a real valued 
function 

Thus M(o) may apply a general modification to yield a 
10 phase equivalent function having a different amplitude to the 
standard correlation function. 

d) M{co) = W x (<u)W*(a))y/{(D) where y/((o) is a real valued 
function 

Thus M(co) may apply a general frequency-dependent scaling 
15 to yield a phase equivalent function having a different 
amplitude to the standard correlation function. Such a 
scaling is also known as emphasis. 

Alternatively, M(co) may be the function D(a>) which is 
the Fourier transformation of the correlation function D(t): 
20 D(t)= £jV x (t + t')W 2 (V)dt< 

D(ca) is mathematically equivalent to W^{co)W*{co) and may 
be arrived at without calculating W x (a>) and W 2 (co)* . This is an 
alternative method to calculating the standard dispersion 
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corrected correlation function. The steps are calculate D(t); 
calculate D(co) and apply a frequency stretching operation to 
arrive at the dispersion corrected correlation function: 

G(t)=— Cb\f(co)]vtv(icQt)d(D . 

5 One advantage of using the dispersion corrected 

correlation function is that it is applicable in situations 
where the precise time, t 0 , at which a contact occurred is not 
known. This is because an offset t 0 (i.e. t 0 *0) in the response 
functions is represented as an additional factor exp(icot 0 ) in 
10 the Fourier transformations, W } (co) and W 2 {cd) which cancels in 

the intermediate function W x (co)W 2 * {of) . 

A transducer may act as both the first and second sensor 
whereby, the dispersion corrected correlation function is an 
autocorrelation function. The autocorrelation function may be 

15 calculated applying the same steps for the dispersion 
corrected correlation function using W x {t)= W 2 (t) . 

The frequency stretching operation may be determined as 
described in relation to Figures la to Id by £l(co) = W[f(co)] in 
which f(co) is chosen so that Cl{co) represents an imaginary 

2 0 response of a medium where bending waves of all frequencies 
travel with identical phase velocity. For example, for bending 
mode vibrations that follow a square-root relation of 
wavevector, k, to angular frequency, k = C^co , the frequency 
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stretching operation f{co) = vC4cj gives an imaginary response, 
Q.(t) , of a medium with constant arbitrary phase velocity v . C 
is a constant and is defined by C=(m/B)^1/4 where /x=mass per 
unit area, and B=bending stiffness. 
5 The inverse Fourier transform, Cl(t) may be examined to 

determine the distance to the site of the contact. Q(Y) , is 
centred at a value t x that is proportional to the distance 
between the site of the contact and each sensor. 

The dispersion corrected correlation function may also be used 

10 to determine the precise location of the contact. For 
example, for the function, G(t) , the centre of the correlation 
function is located at t = {x x -x 2 )/v where x x and x 2 are the 
distances from the contact point to the respective vibration 
sensors. The dispersion corrected correlation function may 

15 therefore be used to determine a difference in path- length 
(i.e. x x -x 2 ) between the contact site and the two sensors. 
This quantity, x x - x 2 , defines a hyperbola of possible 
locations of the contact on the panel surface. If a third 
sensor is provided, a second dispersion corrected correlation 

2 0 function, from a different combination of sensor positions, 
provides a second path-length difference, eg. x x - jc 3 . 
Alternatively, a second pair of sensors may be mounted to the 
member to provide a second dispersion corrected correlation 
function. In many cases (for example when the sensors are 

2 5 positioned in the four corners of a rectangular panel) , the 
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two hyperbolic curves defined by the two path- length 
differences have a unique intersection point which determines 
unambiguously the location of the contact. 

Another advantage of using the dispersion corrected 
5 correlation function is in the treatment of waves reflected 
from boundaries of the member. The reflected waves create the 
effect of virtual sources which are located outside the 
physical boundaries, specifically, at the contact site 
reflected in the axes of the boundaries. The impulse generated 

10 by a contact may show features corresponding to 'echoes' which 
arrive after the direct waves of .the impulse. By applying the 
dispersion correction to an impulse at t=0, discrete 
reflections may be resolved as peaks in the corrected impulse 
response. The location of these peaks may be used as 

15 additional information for determining the location of the 
contact. This may be particularly effective if the 

sensitivity, or acceptance, of the sensors is not independent 
of the direction of the incoming wave since the measured 
signal corresponding to the direct path from the contact may 

2 0 be weaker than the signal due to reflected waves. 

The dispersion-corrected autocorrelation function may be 
applied to locate the path- length difference between the 
sensor and real and virtual sources of bending waves arising 
from reflections at the boundaries. Such information may help 

2 5 the determination of the contact location. 

Reflected waves may be suppressed by placing an absorbing 
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material in contact with the edges of the member. 
Alternatively, the processor may be adapted to remove the 
contribution of reflected waves from the measured bending wave 
signal. This may be achieved on the basis that in a dispersive 
5 medium, i.e. one with a dispersion relation of the form 

k = c4a> , low-frequency components travelling along the direct 
path may arrive after the first high-frequency reflected 
waves. This processor may be adapted to provide a moving 
average low-pass filtering operation on the original measured 

10 bending wave signal where the width of the averaging window 
varies locally with the time coordinate as Atoct 2 . 

The measured bending wave signal may be transformed to a 
response in which the signals due to the direct and reflected 
wave both oscillate with constant but different periods. For 

15 example, a measured bending wave signal W(t) , from a contact 
at time t = 0 may be transformed on the time axis by using 
U(t) = W(l/r) . In the function U{r) the signal from a sharp 
impulse propagating on a dispersive member oscillates with a 
constant period. Furthermore, the oscillations from the direct 

20 wave have a longer period than those of the reflected waves. 
The contribution of reflected waves can therefore be removed 
by a low pass filter operation on U(t) . The response may then 
be transformed back to linear units of time, as W'{t) = U(l/t) . 
The procedure may be generalised to other forms of the 

25 dispersion relation. 
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For a perfect rejection of reflected waves, it is 
necessary to know the dispersion relation in the member, the 
time t 0 at which the contact occurred, and the distance x 
between the contact site and the sensor. However, generally, 
5 only the first is known. Thus the processor may be adapted to 
provide estimates or substitutes for t 0 and x which may be 
substituted into a calculation to remove the reflected waves. 
For example, an estimate for t 0 may be the time at which the 
contact was first detected by any sensor on the member, i.e. 

10 the time at which the measured signal first passes a 
predetermined threshold. The distance x may be set as the 
distance between the relevant sensor and the furthest point on 
the member or the maximum dimension (e.g. diameter) of the 
member. Use of these substitutes should leave the direct -path 

15 signal intact. Estimates of t 0 and x obtained by other means 
may also be used. 

There may be multiple (i.e. n) sensors on the member and 
thus the number of distinct correlation functions is n(n-\)/2 . 
The processor may be adapted to create a mapping function 

2 0 which maps the surface of the member for each correlation 
function whereby the dispersion corrected correlation 
function, G(t), is swept over the surface such that all co- 
ordinates of a given path- length difference, Ax , take the 
value G(Ax/v) . The product of the entire set of mapped 

2 5 correlation functions may then be calculated and the location 
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of the contact may be determined from the maximum co-ordinate. 
This method has the desired property of increasingly rejecting 
spurious contributions from reflections as the number of 
sensors is increased. 

Each sensor and its associated circuitry may have 
identical or different phase properties. Differences in the 
phase properties may be caused by mechanical inconsistencies 
in the manufacture and mounting of the sensors, or by finite 
tolerances of component values in the associated amplification 
and filtering circuits. If the phase properties of each 
sensor differ, the dispersion corrected correlation function 
may be convolved with some asymmetric function. This 
asymmetric function may lead to incorrect measurements of 
path-difference difference since the position of peaks in the 
dispersion corrected correlation function is no longer 
independent of the phase properties of the sensor. 

The processor may thus be adapted to perform the 
following steps : 

a) estimate a convolution correction coefficient (f> 12 (co) from 



where \v hj (G))} and ^V\j{co)} are the Fourier transformation and 
complex conjugate Fourier transformation of two measured 
bending wave signals §V XJ (t)) and {W 2J (t)\ and {Axj\ is the path- 
length difference; 

b) calculate the dispersion corrected correlation function 
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with phase correction from: 

Git) = £ £V, [f (to)] <t> n \f{a>)} exp(iot)da> . 

The phase correction has the effect of deconvolving the 
dispersion corrected correlation function such that the 
position of the maximum is consistent with the position of the 
5 contact . 

The processor may further be adapted to include in the 
determination procedure any available information about where 
the contact can be expected. This may be particularly useful 
in situations where the position of a contact may not be 

10 unambiguously determined by the dispersion corrected 
correlation functions, e.g. when reflected waves interfere 
with the direct -wave signature or less than three sensors are 
used. For example, if the member is an input device for a 
graphical user interface where the user is presented with a 

15 choice of 'buttons' to press, it may be useful to assume that 
any contact on the member occurs within the discrete areas 
corresponding to the buttons . 

Alternatively, a map of the probability at which a 
contact is likely to occur and which is based on the expected 

2 0 behaviour of the user may be used. The device may comprise a 
software application with a graphical user interface (GUI) 
which interacts with the operating system by means of an 
application program interface (API) in which the API is 
adapted to generate the probability map. The probability map 

Atty. Dkt. No. 081-0002-1 
081-0002-1 final patent app.doc 



17 

may be based on the location, size, and frequency of use of 
objects presented by the graphical user interface. The 
probability map may also be based on information about the 
relative likelihood of the various GUI elements being 
5 activated . 

The information in the mapped correlation functions may 
be combined with the probability map to give higher 
reliability. The expected probability may also be a further 
input to a neural -net which accepts information derived from 

10 the sensors. 

The member may comprise a raised pattern on its surface 
whereby a contact drawn across the surface provides a variable 
force to the member to generate bending waves in the member. 
The pattern may be periodic, or quasi -periodic with a 

15 statistically well-defined spatial distribution of 
undulations. The processor may be adapted to determine the 
contact position by using knowledge of the periodicity of the 
pattern and the fact that the interval between impulses 
represents the time in which a contact, which may be provided 

2 0 by a stylus, has travelled to an adjacent feature of the 
pattern. The pattern may be random whereby a contact 

travelling over the surface of the member generates a random 
bending wave signal . 

The use of random surface pattern may be used 

25 independently of the dispersion correction correlation 
function. Thus according to another aspect of the invention, 
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there is provided a contact sensitive device comprising a 
member capable of supporting bending waves, a sensor mounted 
on the member for measuring bending wave vibration in the 
member to determine a first measured bending wave signal and a 
5 processor which calculates information relating to the contact 
from the measured bending wave signal from the sensor, wherein 
a surface of the member comprises a raised pattern whereby a 
contact drawn across the surface provides a variable force to 
the member to generate bending waves in the member. 

10 The device may comprise sensing means to determine a 

second measured bending wave signal which is measured 
simultaneously with the first measured bending wave signal and 
the processor may calculate information relating to the 
contact from a dispersion corrected correlation function of 

15 the two measured bending wave signals. The dispersion 
corrected correlation function is described above and thus 
features of the first and second embodiments may be applied to 
this embodiment. The dispersion corrected correlation 

function is particularly useful if the measured bending wave 

2 0 signals have a large bandwidth centred on a frequency which 
gives a phase-velocity of bending waves in the member which is 
much greater than the maximum lateral velocity of the contact . 

According to another embodiment of the invention, there 
is provided a contact sensitive device comprising a member 

2 5 capable of supporting bending waves, a sensor mounted on the 
member for measuring bending wave vibration in the member to 
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determine a first measured bending wave signal and a processor 
which calculates information relating to a contact from the 
measured bending wave signal from the sensor, wherein the 
device comprises at least a second sensor to determine a 
5 second measured bending wave signal which is measured 
simultaneously with the first measured bending wave signal and 
the processor optimises a product of a set of corrected 
impulse response measurements from each sensor to determine 
information related to the contact . 
10 The information calculated may be the time t 0 at which 

contact occurs and thus the form of the dispersion relation, 
k{co) , and the distance, x, separating each sensor from the 
contact site must be known. The corrected impulse response 
measurement may be calculated by using the following steps: 

15 1) Calculate the Fourier transform W{co) of a measured bending 
wave signal W{t) ; 

2) Calculate an equivalent response, W'{t) , from a notional 
sensor positioned at the contact site by 

W\ai) = W(a))cxp[ik(o))x] . 

3) Calculate the inverse Fourier transform of W f (o)) to provide 

2 0 function JV f (t) . 

The product is thus Y[Wj(t) in which the function W'(t) 

shows an initial impulse from the contact which is without 
dispersion and which is centred at t 0 , the time of the 
impulse. Features in W(i) due to edge reflections will appear 
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later in W'(t) but will not be corrected for dispersion in the 
same way as for the direct impulse which is thus more easy to 
identify. 

The product shows a strong peak due to the direct 
5 impulse, and a weak background noise. By taking the product we 
reinforce the signal coincident at t 0 in all Wj(t) , whereas the 
information due to reflections is uncorrelated and suppressed. 
This procedure is therefore a reliable way of determining t 0 . 

Conversely, the information relating to the contact may 
10 be the location of the contact, i.e. the co-ordinate r and 
thus the sensor-contact distances, Xj and the time t Q for 

which the maximum value of Ylffiffo) is obtained must be known. 

This optimisation process may involve iterative 
refinement of estimates for r and t 0 . The initial estimate may 

15 be derived from impulse response functions whose high temporal 
frequency components have been suppressed, with the result 
that the estimate is easy to obtain (fewer local maxima) , but 
of low spatial precision. Further iterations may introduce 
progressively higher frequency components as the estimate is 

20 refined. 

The various embodiments of the invention therefore 
provides two complementary methods of determining the contact 
position: the method of dispersion-corrected correlation 
functions, and the method of maximising Y1j W j(*o) - 
25 The following characteristics may apply to all 
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embodiments of the invention. The device may comprise means 
for recording measured bending wave signals from the or each 
sensor over time as the contact moves across the member. The 
measured bending wave signals may be recorded as time-series 
5 data, i.e. a sequence of values measured at different times. 

The time-series data may be analysed as a sequence of 
short sections or 'frames' of data, which may have widths or 
durations of 10ms. The processor may apply a window function 
to the frames of data. Window functions are well known in the 

10 art - see for example Press W.H. et al . , Numerical Recipes in 
C, 2 Ed., Cambridge University Press 1992, Ch. 13.4. The 
processor may be adapted to extract information on the contact 
which has been averaged over the duration of the frame, e.g. 
the mean positional co-ordinate of the moving contact. The 

15 processor may apply the dispersion corrected correlation 
technique to calculate the mean path length difference for 
each frame of data to give a sequence of path- length 
differences over time. 

Alternatively, the time-series data may be analysed using 

2 0 adaptive filters such as those described in Grant PM et al 
"Analogue and Digital Signal Processing and Coding", Ch 10 
(1989) . The adaptive filter may correct for a convolution 
function which relates the time-series data from the first and 
second measured bending wave signals. The convolution function 

25 changes slowly over time as the contact moves and is dependent 
on the position of the contact. 
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The adaptive filter may calculate a convolution 
correction which when applied to the measured bending wave 
signals produces a signal as close as possible to the 
measurement input. The first measured bending wave signals 
5 form a first time-series data which may be sent to the 
reference signal of the adaptive filter and the second 
measured bending wave signals form a second time-series data 
which may be delayed, preferably by the maximum expected width 
of the correlation function, before being sent to the primary 

10 input of the adaptive filter whereby the convolution 
correction may be calculated. The processor may apply the 
dispersion correction to the convolution function to give a 
function whose maximum is determined by the difference between 
the path- length from contact to the sensor and the path- length 

15 from the contact to the sensing means. 

The information calculated may be the location of the 
contact or may be other information, e.g. pressure or size of 
the contact. The information relating to the contact may be 
calculated in a central processor. The sensors may be mounted 

2 0 at or spaced from an edge of the member. The sensors may be 
in the form of sensing transducers which may convert bending 
wave vibration into an analogue input signal. 

The member may be in the form of a plate or panel . The 
member may be transparent or alternatively non- transparent , 

2 5 for example having a printed pattern. The member may have 
uniform thickness. Alternatively, the member may have a more 
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complex shape, for example a curved surface and/or variable 
thickness. The processor may be adapted for complex shaped 
members by providing an adaptive algorithm such as a neural 
net to decipher the contact location from the bending wave 
5 signal received by the sensor. 

The device may be a purely passive sensor with the 
bending wave vibration and hence the measured bending wave 
signals being generated by an initial impact or by frictional 
movement of the contact. Alternatively, the device may be an 

10 active sensor and thus the device may comprise an emitting 
transducer. The transducer may have dual functionality, namely 
acting as an emitting transducer and a sensor. The member may 
also be an acoustic radiator and bending wave vibration in the 
member may be used to generate an acoustic output . 

15 Measurements of the bending wave signal due to a contact 

may be contaminated by the bending waves due to the audio 
signal, particularly when the audio signal is similar to the 
bending wave signal generated by the contact . The effect may 
be minimised by ensuring the frequency band of the audio 

2 0 signal differs from and does not overlap the frequency band of 
the measurements from the sensors and sensing means. The 
audio and measured signals may thus be filtered, for example, 
the audio band may be limited to frequencies below 2 0kHz, and 
the vibration measurements may be limited to frequencies above 

25 20kHz. 

The device may be a dual active and passive sensor and 
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may be adapted to switch between active and passive sensing 
modes depending on whether contact is applied to the device. 
The device may cycle between resting in passive sensing mode 
when no contact is detected, switching to active mode sensing 
5 when a contact is applied and returning to passive sensing 
mode once the contact is removed to wait for further contacts. 
This may be advantageous to avoid the power drain when the 
device is in active mode. 

The contact may be in the form of a touch from a stylus 

10 which may be in the form of a hand-held pen. The movement of a 
stylus on the member may generate a continuous signal which is 
affected by the location, pressure and speed of the stylus on 
the member. The stylus may have a flexible tip, e.g. of 
rubber, which generates bending waves in the member by 

15 applying a variable force thereto. The variable force may be 
provided by tip which alternatively adheres to or slips across 
a surface of the member. As the tip moves across of the 
member a tensile force may be created which at a certain 
threshold, causes any adhesion between the tip and the member 

2 0 to break, thus allowing the tip to slip across the surface. 

Alternatively, the contact may be in the form of a touch 
from a finger which may generate bending waves in the member 
which may be detected by passive and/or active sensing. The 
bending waves may have frequency components in the ultrasonic 

25 region (>20 kHz) . Passive sensing is therefore sensitive to 
contacts with both fingers and styli. 
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When the device is acting as an active sensor, i.e. with 
an emitting transducer generating an excitation signal, the 
contact may exert a non- linear force on the member so as to 
generate harmonics of the excitation signal. The processor may 
5 comprise signal processing devices to isolate the excitation 
signal from the harmonics so that the harmonics may used to 
determine the contact position in a similar manner to passive 
sensing. The harmonics effectively constitute a source of 
bending waves from the contact site. 

10 The or each emitting transducer or sensor may be a bender 

transducer which is bonded directly to the member, for example 
a piezoelectric transducer. Alternatively, the or each 
emitting transducer or sensor may be an inertial transducer 
which is coupled to the member at a single point. The inertial 

15 transducer may be either elect rodynamic or piezoelectric. It 
may be possible to use audio transducers which are already in 
place as sensing and/or emitting transducers. 

A contact sensitive device according to the invention may 
be included in a mobile phone, a laptop or a personal data 

20 assistant. For example, the keypad conventionally fitted to a 
mobile phone may be replaced by a continuous moulding which is 
touch sensitive according to the present invention. In a 
laptop, the touchpad which functions as a mouse controller may 
be replaced by a continuous moulding which is a contact 

25 sensitive device according to the invention. The moulding may 
be implemented as a mouse controller or other alternatives, 
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e.g. a keyboard. Alternatively, the contact sensitive device 
may be a display screen, e.g. a liquid crystal display screen 
comprising liquid crystals which may be used to excite or 
sense bending waves. The display screen may present 

5 information relating to the contact. 
BRIEF DESCRIPTION OF DRAWINGS 

The invention is diagrammatically illustrated, by way of 
example, in the accompanying drawings, in which :- 

Figures la to Id are a graphic illustration of a method 
10 of dispersion correction according to the prior art, in which 
Figure la is a graph of a dispersive impulse response showing 
response in arbitrary units against time. Figure lb is a graph 
of a dispersive frequency response showing response in 
arbitrary units against frequency. Figure lc is a graph of a 
15 non- dispersive frequency response showing response in 
arbitrary units against frequency. Figure Id is a graph of a 
non -dispersive impulse response showing response in arbitrary 
units against time; 

Figure 2 is a plan view of a contact sensitive device 
2 0 according to the present invention; 

Figure 3 is a perspective view of a first device 
incorporating passive touch sensing; 

Figure 4 is a perspective view of a second device 
incorporating passive touch sensing; 
25 Figure 5 is a block diagram of a processing algorithm for 

the passive sensing of Figures 3 and 4; 
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Figure 6 is a perspective view of first device 
incorporating active touch sensing; 

Figure 7 is a perspective view of a second device 
incorporating active touch sensing; 

Figure 8 is a flow chart showing a method of calculating 
the path-length difference between the contact and two 
measurement points using the dispersion corrected correlation 
function; 

Figure 8a is a schematic plan view of a device to which 
the method of Figure 8 is applied; 

Figure 8b is a graph of dispersion corrected correlation 
function against time; 

Figure 9 is a flow chart showing a first method of 
removing reflections from a measured bending wave signal, 

Figure 10 is a flow chart showing a second method of 
removing reflections from a measured bending wave signal. 

Figure 11a and lib are schematic perspective and plan 
views of a touch sensitive device according to another aspect 
of the invent ion ; 

Figure 12 is a schematic block diagram of a processing 
algorithm which may be used for the device of Figure 11. 

Figure 12a is a graph showing the combined transfer 
function H(f) against frequency (f) for the filter and 
amplifier of Figure 12; 

Figure 12b is a graph of a measured bending wave signal 
against time; 
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Figure 13 is a flow chart of the steps for obtaining an 
empirical phase correction; 

Figure 13a is a plan view of a grid defined on a member 
for use in the method of Figure 13; 
5 Figure 14 is a touch sensitive device according to 

another aspect of the invention; 

Figure 15 is a circuit diagram of an adaptive noise 
canceller which may be used in the various devices; 

Figure 16a is a schematic block diagram of a contact 
10 sensitive device which also operates as a loudspeaker; 

Figure 16b is a method of separating audio signal and 
measured bending wave signal in the device of Figure 16a; 

Figure 17 is a flow chart showing a method of calculating 
the contact location using the dispersion corrected auto- 
15 correlation function; 

Figure 17a is a schematic plan view of a device to which 
the method of Figure 17 is applied; 

Figure 17b is a graph of dispersion corrected auto- 
correlation function against time, and 
20 Figure 18 is a block diagram showing how an adaptive 

filter may be used to calculate information relating to the 
contact . 

Figure 2 shows a contact sensitive device (10) comprising 
a transparent touch sensitive plate (12) mounted in front of a 
2 5 display device (14) . The display device (14) may be in the 
form of a television, a computer screen or other visual 
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display device. A stylus (18) in the form of a pen is used 
for writing text (20) or other matter on the touch sensitive 
plate (12) . 

The transparent touch sensitive plate (12) is a member, 
5 e.g. an acoustic device, capable of supporting bending wave 
vibration. Three transducers (16) are mounted on the plate 
(12) . At least two of the transducers (16) act as sensors or 
sensing means and are thus sensitive to and monitor bending 
wave vibration in the plate. The third transducer (16) may 

10 also be a sensing transducer so that the system corresponds to 
the passive contact sensitive device of Figure 3 or Figure 4. 

Alternatively, the third transducer may be an emitting 
transducer for exciting bending wave vibration in the plate so 
that the system corresponds to the active sensor of Figure 5. 

15 In the Figure 6 or Figure 7 embodiment, the active sensor may 
act as a combined loudspeaker and contact sensitive device. 

Figures 3 and 4 are more detailed illustration of two 
contact sensitive devices (32,33). The contact sensitive 
devices (32,33) comprises a member in the form of a panel (24) 

2 0 capable of supporting bending wave vibration and three sensors 
in the form of sensing transducers (26) for sensing bending 
wave vibration at their respective mounting points. The 
vibration pattern (28) is created when pressure is applied at 
a contact point (30) . The devices may be considered to be 

2 5 passive contact sensitive devices since the devices do not 
comprise an emitting transducer. Thus the bending wave panel 
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vibration in the panel is generated solely by the contact. 

In a passive sensor an impulse in the body of the panel 
(24) starts a bending wave travelling towards the edge of the 
panel (24) . The bending wave is detected by the three sensing 
5 transducers (2 6) mounted equidistant ly around the edges as in 
Figure 3 or by the three sensing transducer mounted on a 
surface of the panel (24) but spaced from the edges of the 
panel (24) as in Figure 4. The measured bending wave signals 
are processed to determine the spatial origin and force 

10 profile of the applied impulse. 

Figure 5 shows a possible implementation for the 
processing of the bending wave information sensed at each 
sensing transducer (2 6) of Figure 3 or Figure 4. In Figure 5, 
the bending waves in the panel are sensed by three sensing 

15 transducers (26) . The sensing transducers (26) measure 

analogue bending wave signals Wi(t), W 2 (t) and W 3 (t) which are 
transmitted to a multiplexing analogue to digital converter 
(ADC) (54) . The resultant digital input signal is transmitted 
to the central processor (34) from which information (58) 

2 0 relating to the location and profile of the contact impulse is 
determined. 

Figures 6 and 7 are more detailed illustrations of 
alternative combined touch sensitive and audio devices 
(35,37). The devices each comprise a panel (24) capable of 
2 5 supporting bending wave vibration and an emitting transducer 
(31) for exciting bending wave vibration in the panel (24) . 
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The device (35) in Figure 6 further comprises two sensing 
transducers (26) for sensing bending wave vibration at their 
respective mounting points whereas the device (37) in Figure 7 
comprises only one sensing transducer (2 6) . The vibration 
5 pattern (2 8) is interrupted when pressure is applied at a 
contact point (30) . The devices may be considered to be active 
contact sensitive devices since the devices comprise an 
emitting transducer (31) . 

In Figure 6, the sensing and emitting transducers (26,31) 
10 are spaced equidistantly around the edges of the panel (24) 
whereas in Figure 7, the sensing and emitting transducers 
(26,31) are distanced from the edges of the panel (24) and are 
mounted to a surface thereof . 

Figure 8a shows an embodiment having two sensors (102) 
15 mounted on a member (100) to which a contact is applied at a 
contact location (104) . Figure 8 shows a method of calculating 
the dispersion corrected correlation function to reveal the 
difference in path length between the contact location (104) 
and the sensors (102) . The method comprises the following 
20 steps: 

(a) Measure two bending wave signals W x (t) and W 2 (t) ; 

(b) Remove reflections from the measured signals and calculate 
W{(t) and W 2 (t) , e.g. by using the method set out in Figure 9; 

(c) Calculate the Fourier transform of W{{t) and W 2 (t) to 
2 5 arrive at W x (co) and W 2 {cd) and hence the intermediate function 
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W x {co) W*(o)) ; where W*(cd) is the complex conjugate Fourier 
transform. 

(d) and (e) at the same time as performing steps (a) to (c) , 
the frequency stretching operation f(co)=vC^fw is calculated 
5 using the predetermined panel dispersion relation k = C-[a> . 

(f) W } (g>) and W 2 {co) and /{&>)= vC^Jco are combined to arrive at 
the dispersion corrected correlation function: 

Gif) = [jV x [f(co))V 2 * [f(a>)]eKp(ia*)da> ; and 

(g) the dispersion corrected correlation function is plotted 
10 against time with a peak occurring at time t 12 as shown in 

Figure 8b; 

(h) Ax 12 is calculated from t 12 ; Ax 12 is the path- length 
difference between the path lengths x x and x 2 from the first 
and second sensors to the contact . 

15 (i) Axi 2 is used to calculate the location of the contact. 

Alternatively at step (e) , the dispersion corrected 
correlation function with phase correction <|) 12 set out below 
may be used. The calculation of (j) 12 is explained in Figure 13. 

G{t) = [JV X [f{o>W; \f{co)} <f> n \f{co)] exp(icot)dco 

20 Figure 9 shows a method of removing reflections from an 

impulse measurement where the impulse occurs at t=0. The 
method comprises the following steps: 
i) Measure a bending wave signal W x (t) ; 
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ii) Transform the signal on the time axis by using U(r)= W(}/t) ; 

iii) A low pass filter is applied, for example, U' (x) as shown 
to remove all reflected signals. The constant C from the 
predetermined panel dispersion relation k = C-Ja> is used to 

5 define the width of the convolution function At; 

iv) The response may then be transformed back to linear units 
of time, as W'(t)=U(l/t) . 

Figure 10 shows an alternative method for removing 
reflections from an impulse measurement where the impulse 
10 occurs at t=0 . The method comprises the following steps: 

i) Measure a bending wave signal W x (t) ; 

ii) Estimate the distance x between each sensor and the 
location of the contact; 

iii) Use the estimate x and a predetermined dispersion 
15 relation k = c4a> to define an averaging window N(t, t') 

iv) Apply the averaging window N(t, t f ) to the bending wave 
signal W x {t) to remove the effect of reflections. 
N(t,t') is an example of an averaging window with a 

Gaussian shape. A rectangular window can be implemented to 
2 0 give computational efficiency. Multiple applications of a 

rectangular window may be desirable; a large number of 

repeated applications of a rectangular window will produce a 

similar result to N(t,t') . 

Figures 11a and lib show a contact sensitive device (80) 
2 5 comprising a rectangular member (82) capable of supporting 
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bending waves and four sensors (84) for measuring bending wave 
vibration in the member. The sensors (84) are in the form of 
piezoelectric vibration sensors and are mounted on the 
underside of the member (82) , one at each corner. A foam 
5 mounting (86) is attached to the underside of the member and 
extends substantially around the periphery of the member. The 
foam mounting (86) has adhesive surfaces whereby the member 
may be securely attached to any surface. The foam mounting may 
reduce the reflections from the edge of the member. 

10 Two sets of path length differences, Axi 2 =X!-x 2 and 

Ax 3 4=x 3 -X4 are calculated as described in Figure 8; x ± is the 
distance from each sensor to the contact. As shown in Figure 
lib the hyperbolae (85) defined by the two path-length 
differences are plotted and the location of the contact is the 

15 intersection (87) of two hyperbolae. 

Figure 12 is a schematic diagram illustrating the 
implementation of a processing algorithm in the device of 
Figure 11. The sensors (84) measure analogue bending wave 
signals Wi(t), W 2 (t) and W 3 (t) which are passed through an 

20 amplifier and anti-aliasing (low-pass) filter (88) . The 
amplitude of the combined transfer function H(t) of the 
amplifier and anti-aliasing filter is shown in Figure 12a. The 
filtered signals are converted into digital signals by a 
digitiser (90) and stored in a first-in-first-out buffer 

25 having finite length. The buffer comprises two stores, a pre- 
trigger and a post-trigger store (92,93) for signals measured 
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before and after the detection process is triggered 
respectively. 

The central processor (94) determines information relating 
to the location and profile of a contact of the member by the 
following steps: 

a) The central processor performs a threshold test which is 
illustrated in Figure 12b. The measured bending wave 
signal (96) is compared to a predetermined threshold value 
(98) . When the measured signal passes the threshold 
value, the detection process is triggered. 

b) An array of time- series digital input signals is 
transferred from the buffer to the processor. The signals 
include measurements taken before and after the detection 
process is triggered so that a digital measure of the 
entire waveform of the impulse from the contact is 
reconstructed . 

c) The processor shifts the waveform according to the 
estimate of t 0 so that t 0 is set to zero. 

d) The processor removes the effect of reflections from the 
digitised signal as described above with the estimate of 
t 0 as zero and the estimate of x taken as the diagonal 
length of the member. 

e) The processor applies further processing, in particular 
calculating the dispersion corrected correlation function 
for each diagonally opposed pair of sensors and 
calculating information relating to the contact. 
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The further processing applied by the processor may also 
comprise applying a pre -determined phase correction to the 
dispersion corrected correlation function. This may be 
calculated as shown in Figure 13 . 
5 a) Define a grid - for example, for the embodiment shown in 
Figures 11a and lib which has a member (82) which four 
sensors (84) , this may be done by defining points (71) by 
the vector {r-j}. In this example the grid has 8 rows and 7 
columns so the grid is defined by a set of 56 vectors [ri, 
10 r 2 ... r 56 ] . 

b) A user taps on the first point (71) in the grid defined by 
ri and the first and second bending wave signals {Wi f j (t) } 
and {w 2# j(t)} are measured by the first and second sensor of 
each pair of sensors; 
15 c) Step (b) is repeated until the user has tapped on each 
point in the grid; 

d) Calculate the Fourier transforms of the bending wave 
signals ; 

e) Calculate the empirical phase correction: 



20 <f>n(a)) = 



Z ^j(*>)4 ik(a>)Axj 



where Axj is the difference between the path lengths, xi,j & x 2 ,j 
from the first and second sensors to the contact. The path 
lengths are known from the grid coordinates . 

Figure 14 shows a contact sensitive device (70) 
2 5 comprising a member (72) capable of supporting bending waves 
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and three sensors (64) mounted on the member for measuring 
bending wave vibration in the member. A surface of the member 
(72) comprises a raised pattern (66) which is a periodic 
pattern of raised crossed lines. A stylus (78) is drawn across 
5 the surface along a path (74) and as it crosses a line of the 
pattern it generates bending waves (76) in the member. 

Figure 15 shows an adaptive noise canceller for example 
as described in "Widrow et al Adaptive Noise Cancelling: 
Principles and Applications, Proceedings of the IEEE, Vol 63 

10 No 12 pp 1692 (1975)". The adaptive noise canceller comprises 
an adaptive filter (40) which takes the initial audio signal 
as the reference input. The canceller may be used to remove 
the contribution of the audio signal from the output of a 
vibration sensor before any further processing occurs. The 

15 adaptive filter shown is one example of an adaptive filter 
that can be applied to this task. 

Figure 16a shows a contact sensitive device which also 
operates as a loudspeaker. Figure 16b shows a method for 
partitioning the audio signal and measured signal into two 

2 0 distinct frequency bands so that the contribution of the audio 
signal to the processed measured signal is suppressed. The 
device comprises a member (106) in which bending waves are 
generated by an emitting transducer or actuator (108) and the 
contact. The emitting transducer applies an audio signal to 

25 the member (106) to generate an acoustic output. Before being 
applied to the member, the audio signal is filtered by a low 
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pass filter (112) which, as shown in Figure 16b, removes the 
audio signal above a threshold frequency f 0 . 

As shown in Figure 16b, the contact generates a signal 
which has a power output which is substantially constant over 
5 a large frequency band. The signal from the contact and the 
audio signal sum to give a combined signal which is passed 
through a high pass filter (114) to remove the signal above 
the threshold frequency f 0 . The filtered signal is then 
passed to a digitiser (116) and onto a processor (118) . 

10 Figure 17a shows an embodiment having a single sensor 

(120) mounted on a member (100) to which a contact is applied 
at a contact location (104) . Bending waves are reflected from 
the edge of the member and create an image of a virtual source 
which is at location (122) . Figure 17 shows a method of 

15 calculating the dispersion corrected auto-correlation function 
to reveal the contact location (104) . The method comprises 
the following steps: 

(a) Measure one bending wave signal W x (t) ; 

(b) Calculate the Fourier transform of W[{t) to arrive at 
20 W x {(6) ; 

(c) at the same time as performing steps (a) and (b) , the 
frequency stretching operation f{co) = \>c4cQ is calculated using 
the predetermined panel dispersion relation k = c4a> . 

(d) W x (co) and f(co)=vC4co are combined to arrive at the 
2 5 dispersion corrected auto-correlation function: 
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G(t) = £ JT t [f(a>)}v; If {a>)]exp{ia)t)da> ; and 

(g) the dispersion corrected correlation function is plotted 
against time with peaks occurring at time t xi and -tn as shown 
in Figure 17b; 

5 (h) Axu is calculated from tn; Ax i2 is the path-length 
difference between the path lengths Xi and x x 1 from the first 
and second sensors to the contact . 

(i) AX12 is used to calculate the location of the contact. 

Figure 18 shows the adaptive filter of Figure 15 may also 
10 be used to calculate the location of the contact from a device 
comprising two sensors. In general, adaptive filters contain a 
finite-impulse-response (FIR) filter. A FIR filter is 
equivalent to a convolution operation, with some convolution 
function <D (t ) . 

15 The signals W 2 (t) measured by the first sensor are sent 

to the adaptive filter (40) and the signals W 2 (t) measured by 
the second sensor are sent to a delay unit (41) . The delay 
unit delays the signals from the second sensor, preferably by 
the maximum expected width of the convolution function. The 

2 0 delayed signals are then sent to the primary input of the 
adaptive filter. The adaptive filter continually updates the 
convolution function so that an estimate, W 2 {t) , of the primary 
signal input, W 2 (t), may be obtained from the reference input, 
Wx(t) . The convolution operation is defined as follows: 
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The location of the contact is calculated in the 
processor by the following steps: 

a) Extract the FIR convolution function from the internal 
memory of the adaptive filter. 

b) Calculate the Fourier transform of the FIR convolution 
function. 

c) Apply the frequency stretching operation f{ay)=vC4co ; 

d) Calculate the inverse Fourier transform to arrive at F(t). 

F(t) is a phase equivalent of the dispersion correction 
correlation function G(t) and thus the Fourier transforms of 
G(t) and F(t) have equal phase but not necessarily the same 
amplitude. Accordingly, the location of any peaks in the time 
domain for F(t) and G(t) are the same and thus the location of 
the contact may be calculated from the peak of F(t) as 
described above for G(t) . 
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